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Abstract 

We investigate the Hamming networks that their nodes are considered as quantum 
harmonic oscillators. The entanglement of the ground state can be used to quantify 
the amount of information each part of a network shares with the rest of the system 
via quantum fluctuations. Therefore, the Schmidt numbers and entanglement entropy 
between two special parts of Hamming network, can be calculated. To this aim, first 
we use the stratification method to rewrite the adjacency matrix of the network in the 
stratification basis. Then the entanglement entropy and Schmidt number for special 
partitions are calculated analytically by using the generalized Schur complement method. 
Also, we calculate the entanglement entropy between two arbitrary subsets (two equal 
subsets have the same number of vertices) in H(2,3) and i7(2,4) numerically, and we 
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give the minimum and maximum values of entanglement entropy in these two Hamming 
network. 
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1 Introduction 

Entanglement plays a crucial role in quantum information processing, including quantum com¬ 
munication [1,2] and quantum computation [3-5]. It is one of the remarkable features that 
distinguishes quantum mechanics from classical mechanics. 

For decades, entanglement has been the focus of much work in the foundations of quantum 
mechanics, being associated particularly with quantum nonseparability and the violation of 
Bells inequalities [6]. Since entanglement has become regarded as such an important resource, 
there is a need for a means of quantifying it. For the case of bipartite entanglement, a re¬ 
cent exhaustive review was written by the Horodecki family [7] and entanglement measures 
have been reviewed in detail by Virmani and Plenio [8]. One of the operational entanglement 
criteria is the Schmidt decomposition [9-11]. The Schmidt decomposition is a very good tool 
to study entanglement of bipartite pure states. The Schmidt number provides an important 
variable to classify entanglement. The entanglement of a partly entangled pure state can be 
naturally parametrized by its entropy of entanglement, dehned as the von Neumann entropy, 
or equivalently as the Shannon entropy of the squares of the Schmidt coefficients [9-11]. The 
situation simplihes if only so called Gaussian states of the harmonic oscillator modes are con¬ 
sidered [12-17]. The importance of gaussian states is two-fold; hrstly, its structural mathemat¬ 
ical description makes them much more amenable than any other continuous variable system 
(Continuous variable systems are those described by canonical conjugated coordinates x and 
p endowed with inhnite dimensional Hilbert spaces). Secondly, its production, manipulation 
and detection with current optical technology can be done with a very high degree of accuracy 
and control. In [18] the authors quantihed the amount of information that a single element of 
a quantum network shares with the rest of the system. They considered a network of quantum 
harmonic oscillators and analyzed its ground state to compute the entropy of entanglement 
that vacuum fluctuations creates between single nodes and the rest of the network by using 
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the Von Neumann entropy. In [19], jafarizadeh et al, quantify the entanglement entropy be¬ 
tween two parts of the network. To this aim,they compute the vacuum state of bosonic modes 
harmonically coupled through the specihc adjacency matrix of a given network. 

In this paper, we want can calculate the Schmidt numbers and entanglement entropy between 
two special parts of Hamming network by generalized Schur complement method. 

In section II, we give some preliminaries such as dehnitions related to association schemes, 
corresponding stratihcation and The Terwilliger algebra. 

In section III, the generalized Schur complement method is used for calculating the Schmidt 
numbers and entanglement entropy. In this method, we will apply the generalized Schur com¬ 
plement method to the potential matrix in the stratihcation basis several times to calculate 
the entanglement entropy between two arbitrary parts in Hamming network. Then ,by using 
this method, the entanglement entropy and Schmidt number for special partitions are calcu¬ 
lated analytically. Also, we calculate the entanglement entropy between two arbitrary subsets 
(two equal subsets have the same number of vertices) in H{2,3) and H{2,4) numerically, and 
we give the minimum and maximum values of entanglement entropy in these two Hamming 
network. 


2 Preliminaries 

In this section we give some preliminaries such as dehnitions related to association schemes, 
corresponding stratihcation and The Terwilliger algebra. 

2.1 The model and Hamiltonian 

We consider nodes as identical quantum oscillators, interacting as dictated by the network 
topology encoded in the Laplacian L. The Laplacian of a network is dehned from the Adjacency 
matrix as = ktdij — , where ki = Aij is the connectivity of node i, i.e., the number 
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of nodes connected to i. The Hamiltonian of the quantum network thus reads: 

ff = 1(P^P + X^(I + 2gL)X) (2-1) 

here I is the N x N identity matrix, g is the coupling strength between connected oscilla¬ 
tors while = {pi,p 2 , ...,Pn) and = (xi, 0 : 2 ,xat) are the operators corresponding to 
the momenta and positions of nodes respectively, satisfying the usual commutation relations: 
[x,p^] = ihl (we set h = 1 in the following) and the matrix V = I + 2gL is the potential 
matrix. Then the ground state of this Hamiltonian is: 

i>(X) = ^ 2jL)X)) (2-2) 

Where the Ag = ^ jg normalization factor for wave function. The elements of 

the potential matrix in terms of entries of adjacency matrix is 

Vij = (1 + 2gKi)6ij - 2gAij 


2.2 Schmidt decomposition and entanglement entropy 

Any bipartite pure state \'iP)ab G -ff = Ha®Hb can be decomposed, by choosing an appropriate 
basis, as 

m 

\^)ab = ® 1^*) (2-3) 

i=l 

where 1 < m < min{dim{HA);dim{HB)}, and ctj > 0 with = 1- Here |aj) (|&i)) form 

a part of an orthonormal basis in Ha {Hb)- The positive numbers ai are called the Schmidt 
coefficients of \'iP)ab and the number m is called the Schmidt rank of \'iP)ab- 

Entropy of entanglement is defined as the von Neumann entropy of either pA or pB'- 


E ^ -Trp4og2PA ^ TrpBlog2PB ^ -y,Xilog2a^ 


( 2 - 4 ) 
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2.3 Association scheme 

First we recall the definition of association schemes. The reader is referred to Ref.[?], for 
further information on association schemes. 

Definition 2.1 (Symmetric association schemes). Let R be a set of vertices, and let 
Ri{i = 0,1,.be nonempty relations on V (i.e., subset oiV xV ). Let the following 
conditions ( 1 ), ( 2 ), (3) and (4) be satisfied. Then, the relations {Ri}o<i<d onVxV satisfying 
the following conditions 

(1) {Ri}o<i<d is a partition oiV xV 

(2) Rq = (o, o) ! o G id 

(3) Ri = R\ ior t] <i < d, where Rj = {f3, a) : {a, f3) E Ri 

(4) For {a,/3) E Rk, the number pL = \y eV ■. {a,y) E Ri and ( 7 ,/S) E Rj\ does not 
depend on {a,/3) but only on i,j and k, define a symmetric association scheme of class d 
on V which is denoted hy Y = (R, {i?i}o<i<d)- Furthermore, if we have ph = p^j for all 
i,j, k = 0,2, ..., d, then Y is called commutative. 

The number v of the vertices, \V\, is called the order of the association scheme and Ri is 
called i-th relation. 

The intersection number pL can be interpreted as the number of vertices which have relation 
i and j with vertices a and fd, respectively provided that {a,/3) E Rk, and it is the same for 
all elements of relation Rk- For all integers i(0 < i < d), set /tj = p^- and note that Hi 7 ^ 0, 
since Ri is non-empty. We refer to Ki as the Tth valency of Y. 

Let Y = (X, {Ri}o<i<d) be a commutative symmetric association scheme of class d, then 
the matrices Aq, Ai, ..., Ad defined by 


(^■1)0 ,(3 


1 if {a, (3) E Ri, 
0 otherwise 


{a,/3EV) 


(2-5) 



Entanglement entropy 


7 


are adjacency matrices of Y such that 


A, A, = ypip 


( 2 - 6 ) 


fc =0 


From (2-6), it is seen that the adjacency matrices Aq, Ai ,form a basis for a commutative 
algebra A known as the Bose-Mesner algebra of Y. This algebra has a second basis Eq, Ed 
(known as primitive idempotents of F) so that 


Eq — —J, EiEj — 6ijEi, ^— I- 


n 


(2-7) 


i=0 


where, J is the all-one matrix in A. Let P and Q be the matrices relating the two bases for A; 


d 


Ai = ^PijEj, 0<j<d, 


1 


i=0 
d 


Ei = QijAj, 0 <j <d. 


^i=0 


Then clearly 


PQ = QP = nP 


( 2 - 8 ) 


(2-9) 


It also follows that 

AEj = P.jEj, (2-10) 

which shows that the Pij (resp. Qij) is the j-th eigenvalue (resp. the j-th dual eigenvalue ) of Ai 
(resp. Ei) and that the columns of Ej are the corresponding eigenvectors. Thus, m* =rank(i7j) 
is the multiplicity of the eigenvalue Pij of Ai (provided that Pij ^ Pkj for k ^ i). We see that 
^0 = and mi =tracei7j = n{Ei)jj (indeed, Ei has only eigenvalues 0 and 1, so 

rank(i7fe) equals the sum of the eigenvalues). 

Clearly, each non-diagonal (symmetric) relation Ri of an association scheme Y = (C, {Ri}o<i<d 
can be thought of as the network (V, Ri) on V, where we will call it the underlying network of 
association scheme Y. In other words, the underlying network T = {V,Ri) of an association 
scheme is an undirected connected network, where the set V and Ri consist of its vertices 
and edges, respectively. Obviously replacing Ri with one of the other relations such as Ri, 
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for i 7 ^ 0,1 will also give us an underlying network F = {V,Ri) (not necessarily a connected 
network) with the same set of vertices but a new set of edges 

2.4 Stratification 

For an underlying network F, let hF = C" (with n = |V^|) be the vector space over C consisting 
of column vectors whose coordinates are indexed by vertex set of F, and whose entries are 
in C. For all P E V, let \P) denotes the element of W with a 1 in the P coordinate and 0 in all 
other coordinates. We observe {|/9)|/9 G V} is an orthonormal basis for hF, but in this basis, 
W is reducible and can be reduced to irreducible subspaces Wi, i = 0,1,d, i.e., 

fF = Wq © Wi ® ... ® W^i (2-11) 

where, d is diameter of the corresponding association scheme. If we dehne Fj(o) = {P E V : 
{o, P) E Ri} for an arbitrary chosen vertex o E V (called reference vertex), then, the vertex 
set V can be written as disjoint union of Fi(Q;), i.e., 

r = L)r.(a). (2-12) 

i=0 

In fact, the relation (2-12) stratihes the network into a disjoint union of strata (associate 
classes) Fj(o). With each stratum Fj(o) one can associate a unit vector 10^) in W (called unit 
vector of i-th. stratum) dehned by 

1*) = ^ E |a). (2-13) 

V a&Viio) 

where, |q;) denotes the eigenket of a-th vertex at the associate class Fi(o) and Ki = |Fj(o)| is 
called the i-th. valency of the network (kj := ph = |{7 : ( 0 , 7 ) G Ri}\ = |Fj(o)|). For 0 < i < d, 
the unit vectors |0j) of Eq.(2-13) form a basis for irreducible submodule of W with maximal 
dimension denoted by hFo. Since {|0i)})hQ becomes a complete orthonormal basis of hFo, we 
often write 

d 

i=0 


(2-14) 
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Let Ai be the adjacency matrix of the underlying network L. From the action of Ai on 
reference state |0o) (|0o) = I®); with o E V as reference vertex), we have 


^i|0o) — X! 

/9eri(o) 

Then by using (2-13) and (2-15), we obtain 


(2-15) 




(2-16) 


3 Entanglement entropy between two parts of a network 

In order to calculate the entanglement entropy between two equal parts in the graph (half 
hrst strata are in one subset and the other strata are in the second subset), we introduce the 
following process: 

First we want to generalize Schur complement method [19]: Suppose we have the following 


matrix which is composed of block matrices. 


Vj 1 Vj 2 0 


\/ 2 , V 22 ^23 


(3-17) 


y 0 V32 V33 J 

Then we do the generalized Schur complement transformation. This transformation is: 


F| 1 Vi 2 0 


1 0 0 


0 0 


V21 V22 V23 0 1 1/231^33 

0 V32 V33 . ,00 1 


V ,^2 V 22 - V23Vfs%2 0 


0 ^33^1/32 1 


(3-18) 

In our work, we will apply the generalized Schur complement method to the potential matrix 
in the stratihcation basis several times. So in the transformed matrix all of the blocks are 


scalar. Therefore the potential matrix is transformed to a 2 x 2 matrix dually. 


Oil ai 2 


afs a22 


(3-19) 
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The wave function in this stage is 


ij{x,y) = Agexp{--{x y) 


Oil ai2 

fll2 0-22 


f A 


X 




by rescaling the variables x and y: 


1/2 

X = ahi X 


1/2 

y = 0.22 y 


(3-20) 


the ground state wave function is transformed to 


ij{x,y) = Agexp{--{x y) 


(i A 

(A 

(7 ij 

U/ 


where 


- 1/2 - 1/2 
7 — ttii 0120-22 


(3-21) 


(3-22) 


So the ground state wave function is 


ij{x,y) = Age A -4 


(3-23) 


Then we can use following identity to calculate the schmidt number of this wave function, 

1 1 -I- 9f 

—:^((^)^ + (?/)^)) +1— 7^A = (3-24) 


^i/2-^y 2(1- 

In order to calculating the entropy, we apply a change of variable as 

2 


1-E = 


= 


12+1 
12-1 


12+1 


So the above identity becomes 


^exp(-^((xf + (yf)) + (i 2 ^ - lAxy) = Yi’^—^4Anix)'i/Jniy) (3-25) 
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and the reduced density matrix is 


the entropy is 


where pn 


2 / u—1 \n 

u+l\v+l> 


P 


2 

Z/ + 1 


E( 


V -I 
Z/+ 1 


)"'|n)(n 


S{p) = -^p„/o^(pn) 

n 


(3-26) 


(3-27) 


ly-l. 


J2pJog{pn) = + (^)%(77Xt) 


'Z/+ 1' 


(3-28) 


and (n) = 


^(P) = ^/op(^) - ^%(^) (3-29) 

By comparing the wave function (3-24) and the identity (3-25) and dehne the scale we 
conclude that 

u = I X 

(z/^ — 1)^'^^ = —7 X 

After some straightforward calculation we obtain 

V = (3-30) 

1 — 7"^ 

By above discussion we conclude that 

n 

where A„ = (4 j)1/2(^)»/2. 
and the entropy is 


S{p) 


Z/+ 1 
2 


log{ 



V — 



(3-31) 
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3.1 Hamming network 

The symmetric product of d-tuples of trivial scheme with adjacency matrices of /„, 

is association scheme with the following adjacency matrices ( generators of its Bose-Mesner 

algebra) 


Aq — Ifi Iji (^ ... $5 

Ai — ^ ) i^Jn In) ® In ® ••• ® In 

permutation 

A= {Jn-In)®{Jn-In)®-®{Jn-In)®In®--®In ( 3 - 32 ) 

permutation T 

i 

Where Jn is n x n matrix with all matrix elements equal to one. This scheme is the well 
known Hamming scheme with intersection number 

i\ 

bi = i, 1 < f < d, 

Ci = {n — l){d — i), 0 < i < d — 1, ( 3 - 33 ) 

where its underlying graph is the cartesian product of d-tuples of cyclic group Zn- We can 
rewrite the matrix as the following and then calculate the singular value decomposition 

of matrix B 


( 


0 


1 1 


1 




0 B 

B Jn—l In—1 


1 

1 

Jn—l In—1 


\ 


1 
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^ 1 0 ^ 


0 On-l 


V 


0 

Vn - 1 
0 

0 


0 

n — 2 0 

0 -1 


0 

0 

0 


0 0 ... -1 


^ 1 0 ^ 


0 On-l 


(3-34) 


J 


Now we want to introduce new basis proportional to this transformation 


| 0 ) = | 0 ) 


| 1 ) = 


1 


\k) = 


\/n - 1 

n—1 


(|l) + |2) + ... + |n-l)) 




/c = 2, 3,..., n — 1 


y/n-I 

The effect of matrix J„ — on the above basis from equation (3-35) will be 


(3-35) 


{Jn - 4)|0) = Vn- 1|1) 

(J„ - 4)|1) = (n - 2)|T) -h y/n - 1|0) 

{Jn-In)\k) = -\k), k = 2,3,...,n-l (3-36) 

So the matrices (Jn — In)-, {Jn — In)o and (Jn — In)+ are 

(Jn - In)- = yjn- 1|T)(0| 

—n—1 ^ ^ 

0„-4)o = (n-2)|l>(l|-i:|fc>(t 

k=2 

{Jn - In)+ = y/n- 1|6)(T| 

Then we introduce the matrices A_ and Aq as 

A^ = yL<s.. ■ In (Jn In)+ ® In ® ® In 

^ In 0 . . . 0 4 0 (Jn - In)- ® h ® ® h 


(3-37) 
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Ao = ^ /n 0 • • • ® 4 ® (^n - 4)o 0 /n ® • • • «) 4 (3-38) 

In this stage we want to construct the hrst stratum of Hamming network. We begin with 

|0o)i = looTT. 0) = |o) 0 |o) ®® |o) 

By using the equation (3-38), we conclude that 


^-i-l0o)i — 0 


H_|0o)i = y/(n- l)d|0i)i 


So we have 


Va Va i 


I02)i = ^ E |0...0^0...0^0...0) 

yC2 i,j=i i j 


\(j)j)i = ^= E |o...o^o...o^o...o^o...o) 


j — 12 

We are interesting in the effect of A_ on the one stratum of hrst strata |0j)i, 

= Vn- l^(j + l){d - j)\(t)j+i)i 


(3-39) 


(3-40) 


from the above equation we conclude that, without the coefficient i/n — 1, this graph is the 
spin 1/2 representation of the angular momentum operator. 

Also we have 

= j(^-2)|0j)i 

A+\(t)j)i = Vn- l-\/(j + l)(d-j)|0j_i)i (3-41) 

The above basis have only the elements 0, 1. For = 2,3, ...,n — 1. We know that, the 
number of all states are rA. If all elements are consisted from 0,1, we have, 2^ states. Now, 
imagine one of the elements of state become o;*. So, a* can take n — 2 value and can be in 
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Cf place, therefore the number of states with one parameter ctj, is (n — 2)Cf. Therefore,we 
have 2*^“^ states of this kind and each state, is repeated by (n — 2)Cf. So, the number of all 
states with one ctj is — 2)Cf. By the same way, let we have two different ctj and aj, 

then the number of states will be 2^~‘^{n — 2^02- Similarity, if we have d different Q;j,then we 
have (n — 2)*^ states, that all of them are singlet. We can extend the term as following 

n‘^ = 2^ + 2‘^-\n-2)Ct + 2'^-\n-2fC^ + ... + 2‘^-^{n-2rci + ... + 2{n-2Y-^Cti + {n-2Y 

(3-42) 

The hrst term in the above is the number of all states contained only 0 and 1, i.e. there is not 
the parameter ai in them. The second term in the extension is the number of all states which 
have only one a^. In fact the mth term in the above, is the number of all states that these 
states are contained m different ttjS. 

The adjacency matrix for these states with m different ctj is similar to the adjacency matrix 
for the |0i)s, but it will be for d' = d — m, i.e., the effect of operators and A_ on the states 
(|(Pj)) is similar to the binary states {\4>j))- The diagonal elements can be extract from following 
equation: 


Ao\^j)ai,a2,...,am = U {n - 2) - m)\ipj) a,,a2,...,arr, j = 0,1,..., d - m (3-43) 

For other strata, we begin with |0 i)a: = ... 0), then apply the oper- 

j 

ator A_ to this state and obtain the state |02)fc, So 


^ j 

1 

\(j)2)k = -J= - E (w^'^ +w^'^)|0---0 1 0...0 1 0...0) 

yd{d — 2)ii<i2 

1 

UAk = -^= = y (y*i+y*3)|0...0 1 0...0 1 0...0 1 0...0) 

y<i(<i-2)(d-3)/2,.<r<,. T T T 
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1 

^ y + 1 0...0 1 0...0 1 0 ... 0 ) 

Id{d-2)...(d-m) 

\j — lyyir ^]—*2 

(3-44) 

It can be shown that for these states 


A_\(t)m)k = Vn- lJm{d -m - l)|0m+i)fc 


(3-45) 


These states are eqnivalent to the angular momentum 2J = d — 2. For the third stratum we 


il^i2 22 

These states are equivalent to the representation of 2 J = d — A angular momentum. The 
highest weight of these kind of basis, is 


) ki ,^2 V5^r? 




+ . .0 1 0...0 1 0...0 1 0 ... 0 ) 


*l<*2<---<*m = l Ol,«2vi«rJ 


(3-47) 


For the network with an arbitrary d, the number of strata is [|] -|- 1. 

For calculating the entanglement entropy, we use the following equations 


7l_|0j)i = Vn-lJ{j + l){d- j)\(pj+i)i 


7lo|0j)i = j(n-2)|0j)] 


7l+|0j)i = Vn-lJ{j + l)(d-j)|0i-i)i 


(3-48) 
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So the Adjacency matrix of each strata is 


A 


0 c\^ 

c\fd' n — 2 
0 c^2{d' - 1) 


0 0 

c^2{d’ - 1) 0 

2{n - 2) c^3{d' - 2) 


0 

0 


(3-49) 


0 0 
^ 0 0 


c^2{d' - 1) (ci'-l)(n-2) cV^ 

0 c\/^ d'{n — 2) y 


Where c = \/n — 1 and d' is a parameter which is related to the number of strata, for the first 
term in equation (3-48) in the binary Hamming, the d' is equal to d and for the second term 
in binary Hamming the parameter d' = d — 2, therefore in each stage this parameter decreases 
consecutively. 


3.1.1 Entanglement entropy in Hamming network between two equal parts: First 
half and second half strata 

In this section, we want to calculate bipartite entanglement between two equal parts of strata in 
Hamming network, after applying the generalized Schur complement method to the potential 
matrices of each strata in Hamming networks, we have a 2 x 2 matrix hnally, which it’s entries 
are 


012 = -2gc\ ^ ^ )( ^ ^ ^ ) = -gc{d' -h 1) = -gVn- l(d' -h 1) 


Oil = X — ak 


(^k-l 


^ ^k—1 


^k-2 


-“ 3 “ 


Where x = 1 + 2gd{n — 1), 

ai =-2g{i-l){n-2), i = l,2,..., 
(jJi = Ag^c^i^d' — i -f 1), i = 1, 2 ,..., 


d' + l 
2 

d'- 1 


(3-50) 
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(I22 — ^ — 


X 


^3 - 

^ _/ _ 1 


a' = -2g{d' - i + l)(n - 2 ), i = 1 , 2 ,..., 


d' + l 


oj[ = Ag^c^i{d' — i + 1), i = 1, 2,..., 


d' -I 


(3-51) 


Then the parameter 7 ^ ( 7 * for ith stratum) for each strata is obtained from (3-22): 


/aiia22 


(3-52) 


3.1.2 Entanglement entropy in Hamming network between two eqnal parts: even 
strata are in first snbset and odd strata are in second subset 

The entanglement entropy and Schmidt number are considered between two equal parts of 
hamming network, that even strata are in hrst subset and odd strata are in second subset.In 
this kind of partitioning we can consider for two case. First case for d = odd and second case 
for d = even. 

Case I((i is odd):In this case the potential matrix Vn can be written as following: 


1 -|- 2gdc^ 


l + 2g{dc^-2{c‘^-1)) ... 


(3-53) 


0 1 -|- 2g{dc^ — {d — l)(c^ — 1 )) 


d-|-l 
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And V 22 is 


1 + 2g{dc^ - - 1)) 0 

0 1 + 2g{dc^ — 3{c^ — 1)) 


0 

0 


V 


0 1 + 2g{dc^ — d{c^ — 1)) J 


And the connection matrix of potential matrix V 12 is: 

^ —2gcVd 0 

-2gc^2{d-l) -2gc^3{d-2) 


0 

0 


Therefore,we have parameters 7 


7i = 


-2gcp{d-l) -2gcVd)^ 




72 = 


+ 2g{n{d - + 2g{n{d - 1 ) - (n - 2 )) 

_ (2 + %V^^ _ 

+ 2g{n{d - 1) - 2(n - 2))^1 + 2g{n{d - 1) - 3(n - 2)) 


{2 + 4mgV^r^ 


7d+i = 

^ A + 2g{n{d — 1) — {d — l){n — 2))Jl + 2g{n{d — 1 ) — d{n — 2 )) 


Where c is 

Case II((i is even):In this case the potential matrix Vn is 




1 + 2gdc^ 0 

0 l + 2g{dc^-2{c^-1)) ... 


0 

0 


0 1 + 2g{dc^ — d{c^ — 1 )) j 


|+lx|+l 


+ 1 w d-\-l 
2^2 

(3-54) 


(3-55) 


(3-56) 


(3-57) 


V 


0 
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And V 22 is 

/ 


1 + 2g(dc^ - (c2 - 1)) 0 

0 1 + 2g{dc^ — 3(c^ — 1)) 


0 

0 


0 1 + 2g(dc^ — (d — l){c^ — 1)) j 


And the connection matrix of potential matrix V 12 is: 


( 


—2gc\fd 0 

-2gc^2{d-\) -2gc^‘i{d-2) 


0 

0 


V 


Therefore, we have | parameters 7 


-2gc^2{d — 1 ) —2gcVd j 


i+ixf 


7i = 


Ag^/rl^^ 


72 = 


+ 2g{n{d - 1))^1 + 2g{n{d - 1) - (n - 2)) 

(4 + 4:)g^/Tl^^ 

1/1 + 2g{n{d - 1) - 2(n - 2))^1 + 2g{n{d - 1) - 3(n - 2)) 


d ^ d 
2^2 

(3-58) 


(3-59) 


{2 + 4{l-l))gV^^ 


7d+i = , . ■ (3-60) 

^ 1/1 + ‘^ 9 {n{d — 1) — (d — 2)(n — 2))y 1 -|- 2g{n{d — 1) — {d — l){n — 2)) 

Finally the entropy of entanglement of each strata, i.e. S{p), is obtained from Eq.(3 — 29). So 
the total entropy is S{p) = EiS{pi). 


3.1.3 Entanglement entropy in Hamming network between two eqnal parts of 
adjacency matrix 

In this part, we want to calculate the entanglement entropy in Hamming network between two 
equal parts of adjacency matrix. 
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Case I(n is even): we can define the blocks of adjacency matrix as following: 

An = A 22 = + (3-61) 

And, the connection matrix is 


Ai 2 = Jll®In® ••• ® In 
So, we have the d types of parameter 7 


(3-62) 


71 

72 

73 


ng 

1 + ng 
ng 

1 -I- “ing 
ng 

1 -|- bng 


Id = 


ng 


(3-63) 


1 -|- {2d — l)ng 

For d = 2, the degeneracy number of 71 is 1 and the degeneracy number of 72 is n'^~^ — 1. 

For d > 3, the degeneracy numbers of 7 * respectively are 

1 , (d — l)(n — 1) , (d — 2)(n — 1)^ , (d — 3)(n — 1)^ , ... , (d — (d — l))(n — 1)'^“^. 

Case II(n is odd):we can dehne the blocks of adjacency matrix as following: 

An = ( J n+l — I n+l ) 0 ® ... 0 ~l~ I n+1 ® {Jn — dn)0...0dn“l“...“l“ In+1 ® In® ■■■ ® In® {Jn — In) 

22 2 2 

(3-64) 

A 22 ~ {J Uni — I n-l ) 0 0 ... 0 Ai -|- I n-1 ® {Jn — In)®---®InI~-"~\~ I n-l ®In®---®In®{Jn — In) 

22 2 ' 2 

(3-65) 

And, the connection matrix is 


Ai 2 — J n+l ,, n-l ® In ® ■■■ ® In (3-66) 

Where Jn+i^n^ is the matrix with rows and columns, with all elements are 1. 

So, we have the d types of parameter 7 

J{n + l){n-l)g 
0^1 “ — ^ — - 
^l + (n + l)^iyi + (n-l)^ 



Entanglement entropy 


22 


_ \/{n + l){n - l)g 

_ \/{n + l){n - l)g 

Jl + (5n + l)gJl + (5n - 1)^ 


For d = 2, the degeneracy number of 71 is 1 and the degeneracy number of 72 is — 1. 

For d > 3, the degeneracy numbers of 7 * respectively are 

1 , (d — l)(n — 1) , (d — 2)(n — 1)^ , (d — 3)(n — 1)^ , ... , (d — (d — l))(n — 1)'^“^. Finally the 
entropy of entanglement of each strata, i.e. S{p), is obtained from Eq.(3 — 29). So the total 
entropy is S{p) = EiS{pi). 


'1 + {{2d-l)n-l)g 


I J 


'1 +{{2d-l)n + l)g 


3.2 Entanglement entropy between all kinds of two parts in Ham¬ 
ming H(2,3) network 

In this section we separated Hamming id(2, 3) network into to parts,the first part contains 5 
nodes and the second part contains 4 nodes, then calculated the entanglement entropy between 
all kinds of bisection numerically. In this case, there are 126 kinds of partitioning in a way that 
two parts. Numerical calculations show that some of these 126 sets have the same entropy, 
such that there are 5 kinds of different values for entanglement entropies. 
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FIG I: Hamming iJ(2,3) graph. 

The maximum entanglement entropy in this section is for subset: 1 of table 1. In this case 
there are the maximum edges between two parts. The minimum entanglement entropy in this 
section are for subset:5 of table 1. In these cases there are the minimum number of edges 
between two parts. 


3.3 Entanglement entropy between all kinds of two parts in Ham¬ 
ming H(2,4) network 

In this section we separated Hamming H{2, 4) network into to equal parts, then calculated 
the entanglement entropy between all kinds of bisection numerically. The order of name of 
vertices are based on H = /4 0 (J 4 — I 4 ,) + — I 4 ) 0 h. 

In this case, there are 6435 kinds of partitioning in a way that two parts. Numerical 
calculations show that some of these 6435 sets have the same entropy, such that there are 22 
kinds of different values for entanglement entropies. The maximum entropy is in the partition 
that the vertices ( 1 , 2 , 5, 7,10,12,15,16) in the hrst part. And, the minimum entropy is in 
the partition that the vertices (1, 2, 3,4, 5, 6 , 7, 8 ) in the first part and other vertices is in the 
second part. In the table 2, the assortment of entropy from maximum to minimum is given 
respect to abundance and agent of each group. 
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set 

partitions 

1 

(2, 3,4, 5, 6), (2, 3,4, 5, 7), (2, 3,4, 5, 8), (2, 3,4,5, 9), (2, 3,4, 6, 9), (2, 3,, 4, 7, 8), (2, 3,5, 6, 9), 

(2, 3, 5, 7, 8), (2,4, 5, 6, 9), (2,4, 5, 7,8), (3,4, 5, 6, 9), (3,4, 5, 7, 8), (3, 5,6, 7, 8), (2,4, 7, 8, 9), 

(3,4, 6, 8, 9), (2, 5, 6, 7, 9), (1, 5, 6, 7,9), (1, 5, 6, 7, 8), (1,4, 7,8, 9), (1,4,6, 8, 9), (1, 3,6, 7, 8), 

(1, 3, 6, 8, 9), (1, 2, 7, 8, 9), (1, 2, 6, 7,9), (1, 2,4, 7, 8), (1, 3, 5, 7, 8), (1, 3,5, 6, 9), (1, 3,5, 6, 7), 

(1,3,4, 8, 9), (1,3,4, 6, 9), (1,3,4, 6,8), (1,2, 5, 6, 7), (1,2, 5,6, 9), (1,2,5, 7, 9), (1,2,4, 7, 9), 

(1,2, 4, 8, 9), 

2 

(2, 3,4, 8, 9), (2, 3, 5, 6, 7), (2,4, 5, 7,9), (3,4, 5, 6, 8), (1, 6, 7,8, 9), (1, 3,4, 7, 8), (1, 3,5, 6, 9), 

(2,3,4, 8,9),(1,2,4, 6,9), 

3 

(3,4, 5, 8, 9), (3,4, 5, 6, 7), (2,4, 5, 8, 9), (2,, 4, 5, 6, 7), (2, 3, 5,6, 8), (2, 3, 5, 7, 9), (2, 3,4, 7, 9), 

(2, 3,4, 6, 8), (2,4, 6, 7, 9), (2,4, 6, 8, 9), (2,, 5, 6, 7, 8), (2, 5, 7,8, 9), (3,4, 6, 7, 8), (3,4, 7, 8, 9), 

(3, 5, 6, 7, 9), (3, 5, 6, 8, 9), (1, 5, 7, 8,9), (1, 5, 6, 8, 9), (1,4, 6, 7, 9), (1,4,6, 7, 8), (1, 3, 7, 8, 9), 

(1,3, 6, 7, 9), (1,2, 6, 7, 8), (1,2, 6, 8,9), (1,2, 3, 7, 8), (1,4, 5,6, 9), (1,4,5, 7, 8), (1,2,3, 6, 9), 

(1,3,4, 5, 6), (1,3,4, 5, 8), (1,2,4, 5, 7), (1,2,4, 5, 9), (1,2, 3,5, 7), (1,2,3, 5, 6), (1,2,3,4, 9), 

(1,2,3,4,8) 

4 

(3,4, 5, 7, 9), (2,4, 5, 6, 8), (2, 3, 5, 8,9), (2, 3,4, 6, 7), (4, 5, 7,8, 9), (4, 5,6, 8, 9), (4, 5,6, 7, 9), 

(4, 5, 6, 7, 8), (2, 3, 7, 8, 9), (2, 3, 6, 8,9), (2, 3, 6, 7, 9), (2, 3, 6, 7, 8), (5, 6, 7, 8, 9), (4, 6, 7, 8, 9), 

(3, 6, 7, 8, 9), (2, 6, 7, 8, 9), (1, 2,4, 6,8), (1,4, 5, 6, 8), (1,4, 5, 7, 9), (1, 3,5, 7, 9), (1, 3,5, 8, 9), 

(1,3,4, 7, 9), (1,3,4, 6, 7), (1,2, 5, 6,8), (1,2, 5, 8, 9), (1,2, 3,6, 7), (1,2,3, 8, 9), (1,2,4, 6, 7), 

(1, 3,4, 5, 7), (1, 3,4, 5, 9), (1, 2,4, 5,6), (1,2,4, 5, 8), (1,2, 3,5, 8), (1,2,3, 5, 9), (1, 2,3,4, 7), 

(1,2,3,4,6) 

5 

(2,4, 6, 7, 8), (2, 5, 6, 8, 9), (3,4, 6, 7,9), (3,4, 7, 8, 9), (1,4, 5,8, 9), (1,4,5, 6, 7), (1, 2,3, 6, 8), 


( 1 , 2 , 3 , 7 , 9 ),( 1 , 2 , 3 , 4 , 5 ) 


Table 1: All equal subsets in H{2,3). The entanglement entropy for these subsets are Si > 
S2>S3>S^> S^. 
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set 

abundance 

agent of partitions 

1 

36 

(1,2, 6, 7,10,12,15,16) 

2 

9 

(1,2,5,6,11,12,15,16) 

3 

288 

(1,2,3,5,8,10,12,15) 

4 

288 

(1,2,3,5,6,11,12,16) 

5 

576 

(1,2,3,5,6,9,12,15) 

6 

576 

(1,2,3,5,6,9,12,16) 

7 

432 

(1,2,3,5,6,8,11,16) 

8 

288 

(1,2,3,5,6,9,11,16) 

9 

24 

(1,2,3,5,6,7,12,16) 

10 

576 

(1,2,3,5,6,8,9,15) 

11 

288 

(1,2,3,5,6,7,9,16) 

12 

288 

(1,2,3,4,5,6,11,16) 

13 

72 

(1,2,3,5,6,8,9,14) 

14 

1152 

(1,2,3,4,5,6,9,15) 

15 

288 

(1,2,3,4,5,6,11,15) 

16 

72 

(1,2,3,4,5,6,11,12) 

17 

288 

(1,2,3,4,5,6,9,14) 

18 

288 

(1,2,3,4,5,6,9,11) 

19 

96 

(1,2, 3,4, 5, 6, 7,12) 

20 

360 

(1,2, 3,4, 5, 6, 7, 9) 

21 

144 

(1,2,3,4,5,6,9,13) 

22 

6 

(1,2, 3,4, 5, 6, 7, 8) 


Table 2: agent of equal subsets in if (2,4) 


















































Entanglement entropy 


26 


4 Conclusion 

The entanglement entropy is obtained between two parts in the Hamming networks that their 
nodes are considered as quantum harmonic oscillators. The generalized Schur complement 
method is used to calculate the Schmidt numbers and entanglement entropy between two 
parts of Hamming graph. Analytically, entanglement entropy in two special partitions are 
calculated and the maximum entropy in different partitioning is given. Numerical results show 
that, the maximum and minimum entanglement entropy have a relation with the edges between 
two parts. 

One expects that the entanglement entropy and Schmidt numbers can be calculated in Johnson 
networks. 


5 Appendix 


5.1 Other expression for stratification of Hamming graph 


In this section,we want obtain the matrix Aij by straight way.We know the basis of different 
strata is 

I00...0), I10...0), I01...0),..., I00...1),..., I11...1) (5-68) 


So,we know basis of each strata from number of l.in m’th strata,the basis are:ej^^i 2 ,...,*m) where 
A < A < ••• < fm-Therefore, the adjacency matrix is dehned as following 

m m 

Jl j'm+l *^*272 " "" JT ( ^ ^ikjm+l) A ^ ^ikjk) ■■■ 

k = l k=l 

(5-69) 

From the equation (3 — 32), The adjacency matrix in Hamming graph commutate with per¬ 
mutation. If TT is permutation, we have 


[A, tt] = 0 


(5-70) 
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It means that,the adjacency matrix under displacement of particles is invariant. 


Let l"^) is the state, that is invariant under permutation of vr^+i 

'Tm+lIV’) = ^|V’> 


(5-71) 


(5-72) 


That 0 is the phase.So, the acting of A on the \ip) is the eigenvector of 7t„j.so, that is enough 
hnding the eigenvectors of permutation group.And, we know that the eigenvectors of permu¬ 
tation group is discrete fourier transform.So,we have 


21 22 


E 1 0 ... 1 0 .. 


a\,a2,...,OLr^ 


y/mld 

Where e is Levi-Civita symbol 


(5-73) 


5.2 Finding the normalization factor for strata of Hamming net¬ 
work 


Let 

21 22 

lS>t = i: (a.;" -ncf)|0,..0'd' 0,..0'd' 0.,.0) (5-74) 

il<*2 

Now,we must hnd the normalization multiplier of 102 )a:-T herefore, we de£ne|02)A:' orthogonal 


to |02)fc- 


Therefore 


h *2 

\h)k' = E +^f^)|o...o'1 "o-.-q'T0...0) 

n<i2 

t(SlS)t. = i E (‘s'?' + = 

^{{d - 2)cn(^-^Ti + = d{d - 2) 


(5-75) 


(5-76) 
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So, the |02)fc is 




1 O-O 1 0-0) 


■\J d{^d 2y ii<i2 

finally, by the same way,we have 


i\ 22 im 

\^)k = E + - + a;"*™)|0...o"E 0...o"E 0...o"E 0...0) 


and 


21 22 


Therefore 


(5-77) 


(5-78) 


\4',n)t'= y (w;‘‘+o;“ + ...+in“”)|0...0 1 0...0 1 0...0 1 0...0) (5-79) 

2l <^22 ‘C • • •'^im 


k{(j>m\(j>m)k' — - 7 X! + — 


So, the \(j)m)k is 

10 m) A: 


khf^k'h , r^k'i 2 , , ^k'im^ _ d{d - 2){d - 3)...{d -m- l){d - m) 


{m — 1)! 


(5-80) 


(m — 1)! 


21 22 


\ d{d — 2)...{d — m) ■ 


E {uj^"^+u^"^ + ...+u’^^^)\ 0...0 1 0...0 1 0...0 1 0 ... 0 ) 


*l<*2<---<*n 


(5-81) 
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